Expansion of the Vortex Cores in YBa2Cu306.95 at Low Magnetic Fields 
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Muon spin rotation (/^SR) spectroscopy has been used to 
measure the effective size ro of the vortex cores in optimally 
doped YBa2Cu3 06.95 as a function of temperature T and 
magnetic field H deep in the superconducting state. While 
ro at = 2 T is close to 20 A and consistent with that 
measured by STM at 6 T, we find a striking increase in ro at 
lower magnetic fields, where it approaches an extraordinarily 
large value of about 100 A. This suggests that the average 
value of the superconducting coherence length in cuprate 
superconductors may be much larger than previously thought 
at low magnetic fields in the vortex state. 
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The vortex core in a type-II superconductor is a region 
in which the superconducting order parameter V'(r) is 
strongly suppressed. The size of the vortex core is there- 
fore closely related to the coherence length ^ — which is 
the smallest length over which -0(r) can change apprecia- 
bly. Some time ago Caroli et al. Q predicted that a dis- 
crete spectrum of quasiparticle excitations existed within 
a radius ^ of the vortex axis. An STM experiment on the 
conventional type-II superconductor NbSe2 by Hess et 
al. j2| confirmed the existence of these localized states in 
the vortex core. Since then both ^SR H and STM |] 
measurements have shown that the vortex core size ro 
in NbSe2 decreases with increasing magnetic field — in a 
manner which scales with the increased strength of the 
vortex-vortex interactions. These same techniques [^^|J^] 
have also shown that tq shrinks with decreasing tempera- 
ture, as predicted by the so-called "Kramer-Pesch effect" 

The electronic structure of the vortex cores in the 
high- Tc superconductors (HTSs) is less certain. Exper- 
iments performed on YBaiCugOy^a (YBCO) |,§ and 
Ndi.85Ceo.i5Cu04_5 (NCCO) seem to support the 
existence of a few bound quasiparticle states in the vor- 
tex core. On the other hand, no evidence of such lo- 
calized states was found in a recent STM study of the 
vortex core in under and overdoped Bi2Sr2CaCu208+5 
(BSSCO) Instead the tunneling spectra in the vor- 
tex core resembled the pseudogap which forms at the 
Fermi surface in the normal state of this material. It 
should be noted that the vortex lattice studied in Ref. |ll| ] 
was either melted or consisted of randomly pinned pan- 
cake vortices. 

It is now widely accepted that the order parameter 
in the hole-doped HTSs (which excludes NCCO) have 
a (j^2_„2-wave symmetry. Several authors have pointed 
out H^jl^ that bound quasiparticle states are unlikely 
to exist in a vortex of a d^2_y2-^&\e superconductor be- 
cause of the nodes which are present in the energy gap 
function A^ = Ao(fc^— fc^) along the directions |fcx| = |fcj/|- 
One way of explaining experiments on HTSs which sup- 
port localized states in the vortex core, is to introduce 
additional components into ■(/;(r) {e.g. a (ia;j,-wave compo- 
nent PI ). However, there is currently no direct evidence 
for an order parameter of mixed symmetry in the bulk of 
HTSs. 

In a previous study p^], we determined the size 
of the vortex cores in the underdoped compound 
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YBa2Cu3O6.60- The electronic structure of the cores 
in underdoped YBCO has yet to be investigated with 
STM — although one may anticipate a local density of 
states resembling the normal state pseudogap, as in un- 
derdoped BSCCO In Rcf. ro was found to 
change as a function of T and iJ in a manner simi- 
lar to that observed in NbSe2, but with a considerably 
weaker temperature dependence. Hayashi et al. ||l6| re- 
cently suggested that this may indicate that the quan- 
tum limit is established at a much higher temperature 
in YBa2Cu3O6.60 because of the small value of ^ relative 
to that in NbSe2. Generally speaking, ^ is considered to 
be "short" in the HTSs {e.g. < 20 A) 0. This is one 
of the primary features which distinguishes them from 
conventional superconductors. In YBa2Cu3O6.60, ''o was 
found to be as large as 80 A at low magnetic fields, which 
suggests that ^ is in excess of 20 A. However, this is a 
phase with a Tc of 60 K, so ^ is expected to be somewhat 
larger than in a higher sample. Thus, it is of great 
interest to check the size and field dependence of rg in 
optimally doped YBa2Cu3 06.95, with the maximum 
of 93 K. Unlike in underdoped YBCO, the formation of 
the pseudogap in YBa2Cu3 06.95 appears to coincide with 
the superconducting gap at Tc. The STM study of the 
vortex core by Maggio-Aprile et al. ^ in near optimally- 
doped YBCO, clearly shows two peaks within a gap like 
structure — indicating the presence of bound quasiparticle 
states. 

In this Letter we present /iSR measurements of the ef- 
fective vortex core size in YBa2Cu3 06.95 as a function of 
temperature and magnetic field applied along the c axis 
of the crystals. Our measurements of tq are shown to be 
consistent with the STM study of Ref. gl at if = 6 T. 
Surprisingly, we find that at low fields ro increases to a 
comparatively large value of 100 A. A simple interpre- 
tation of this result would be that the length scale over 
which the order parameter changes in the region of the 
vortex core {i.e. the definition of the coherence length as 
pertaining to the vortex core) is considerably larger than 
the nominal and accepted value of about 20 A measured 
in high magnetic fields. 

We report here measurements of rg in three different 
samples of YBa2Cu3 06.95. The magnetic field distribu- 
tions in two of these samples were previously recorded 
psf . The first (TWl) was a mosaic of three crystals 
whereas the other (TW2) was a single crystal. All of 
these crystals contained twin boundaries and had tran- 
sition temperatures of 93.2(0.25) K. We also report here 
measurements taken in TW2 after removing the twin 
boundaries. Detwinning was achieved by applying uni- 
axial stress to the crystal with the sample heated to no 
more than 250 °C in an oxygen atmosphere. Subsequent 
to the mechanical detwinning process, the crystal was 
reannealed to set the oxygen doping level. The /j,SR ex- 
periments were performed on the M15 and M20 surface 
beamlines at TRIUMF. Our experimental setup is de- 
scribed elsewhere p9| ]. 

Although /zSR does not directly probe the electronic 



structure of the vortex cores, it does sample the distribu- 
tion of local magnetic fields in the vicinity of the cores. 
The spin of an implanted muon precesses at a frequency 
which is directly proportional to the local magnetic field 
at the muon site. Since the local magnetic field rises to 
a maximum in the vortex cores where superconductivity 
is destroyed, ro is directly related to the high-field tail in 
the measured field distribution. The size of the core is not 
strictly defined however, because there is no sharp discon- 
tinuity in spatial quantities between a normal vortex core 
and the surrounding superconducting material. Here, as 
in our previous work 15| , we define ro to be the radius 
about the vortex axis at which the supercurrent density 
Js (r) reaches its maximum value. This feature allows us 
to accurately monitor changes in the effective size of the 
vortex cores. The supercurrent density Js(r) is obtained 
from the field profile B(r) through the Maxwell relation 
Js(r) — I VxB(r)|. In fitting the measured muon spin pre- 
cession signal, some modelling of B(r) is required. How- 
ever, to appreciate the accuracy of the present study it 
is important to realize that the Js (r) profile does not de- 
pend on the validity of the model assumed, since it is 
essentially the same for any function B(r) which fits the 
data well. 

As in Refs. [^,^, the local field due to the vortex 
lattice at any point in the a-b plane was modelled with a 
theoretical field distribution generated from a Ginzburg- 
Landau (GL) model M 
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Bq is the average magnetic field, G are the reciprocal lat- 
tice vectors, b = Bq/ Bc^, £,ab is the GL coherence length 
and Ki{u) is a modified Bessel function. We do not ex- 
pect the conventional GL model to be valid deep in the 
superconducting state. However, this model gives a very 
good fit to the measured field distribution — which is all 
that is required to generate the corresponding Js (r) pro- 
file needed to determine ro. The summation in Eq. (^ is 
taken over all reciprocal lattice vectors G of a triangular 
vortex lattice. This assumption is reasonable, because for 
field-cooled samples the vortex lattice geometry at low T 
is governed by the geometry of the lattice at the pinning 
temperature. We have shown previously |Q that the 
vortex lattice in our YBa2Cu306.g5 crystals is strongly 
pinned at low T and remains so upon warming up to 
T«0.7 Tc. The lattice in the detwinned crystal depins 
at similar high temperatures. In a dx^-y2-wa?fe supercon- 
ductor, the vortex lattice is predicted to be nearly trian- 
gular at temperatures close to Tc {e.g. see Ref. [^). We 
note that even if the lattice is not triangular, a good fit 
still yields the appropriate Js(r) profile. 

All of the data were fit in the time domain with a the- 
oretical muon polarization function constructed from the 
field profile of Eq. (@). This was multiplied by a Gaussian 
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relaxation function e~'^ * to account for any residual 
disorder in the vortex lattice and the contribution of the 
nuclear dipolar moments to the internal field distribu- 
tion. The residual background signal was fit assuming a 
Gaussian broadened distribution of fields. The Fourier 
transform (FT) of the muon precession signal approxi- 
mates the internal field distribution and resembles the 
predicted asymmetric lineshape for an ordered lattice of 
vortices. However, the FT suffers from noise and broad- 
ening effects associated with the finite number of events 
and limited time range, so that the data must be ana- 
lyzed in the time domain. 

Figure 1(a) shows the first 1.5 /j,s of a typical muon pre- 
cession signal in YBa2Cu3 06.95 displayed in a reference 
frame rotating at 3.3 MHz below the Larmor precession 
frequency of a free muon. The solid curve is a fit (ac- 
tually performed over the first 6 /is) to the theoretical 
polarization function assuming ^afc = 54 A, with all other 
fitting parameters unconstrained. The difference between 
this fit and the measured spectrum is shown in Fig. 1(b), 
compared to the same assuming £^ab = 20 A. Note that 
the quality of the fit is most affected by a change in S^ab 
at early time, where the amplitude of the signal origi- 
nating from the vortex lattice is largest. The ratio of 
to the number of degrees of freedom (NDF) for fits 
assuming different values of ^afc is shown in Fig. 2(a) for 
two different magnetic fields. Due to the high statistics 
of the measured field distribution [i.e. typically consist 
of 2 X 10^ muon decay events), x^/NDF is greater than 
one. Figure 2(b) shows that K = Xab/£.ab obtained from 
the same fits also depends on H. 

Figure 3(a) shows the temperature dependence of tq 
at iJ = 0.5 T and 1.5 T in sample TWl. The error bars 
represent the statistical uncertainty in the fitted values 
of ^ab- Both sets of data show a slight decrease in rg with 
decreasing temperature which is essentially linear below 
50K. The strength of the term linear in T is compara- 
ble to that previously reported in YBa2Cu3O6.60 [ P^ , 
which indicates that thermal vibrations are an unlikely 
source of the observed behaviour. On the other hand, 
this linear term is considerably weaker than in NbSe2 
— a result consistent with the prediction of Ref. ||l^] 
that a vortex core containing only a few discrete bound 
states will stop shrinking below a relatively high satu- 
ration temperature. It should be noted, that the core 
is also predicted to shrink with decreasing T in a pure 
da.2_j^2-wave superconductor ||l^, which may not contain 
any localized quasiparticle states. 

Figure 3(b) shows the magnetic field dependence of 
rg extrapolated to T = 0. Included in this figure is data 
from all three samples. Note that there is excellent agree- 
ment between the twinned and detwinned crystals. In 
Ref. it was shown that detwinning a YBa2Cu3O6.60 
single crystal had virtually no effect on the deduced val- 
ues of To, even though there was evidence for distortions 
in the vortex lattice caused by twin boundary pinning. 
Wc find that the fitted parameter £^ab changes as a func- 
tion of T and H in a manner similar to that of , with 



fa6«[0.828(23)ro-H0.72(1.24)] A. Thus at large H where 
ro is small, we have S.ab ^ tq- It is important to note 
however, that unlike tq, the precise behaviour of ^ab does 
depend on the theoretical model used for B{r). Although 
the GL model is invalid here, our measurements in both 
NbSe2 and YBa2Cu307_5 qualitatively resemble the T 
and H dependence of the core size predicted from a qua- 
siclassical treatment of a vortex in both an s-wave and a 
da;2_j,2-wave superconductor p^ , p3[ |. 

Some of the data in Fig. 3(b) is replotted in Fig. 4, 
along with the measurements of ro{H) in YBa2Cu3O6.60 
(from Ref. [^). As was the case in NbSe2, both Xab 
and the ratio Xab/^o (~k) at r = are well described by 
relations linear in H, so that 



ro{H) = 



Kb{H) 



Xab{H)/roiH) 



Recently the field dependence of Xab has been attributed 
to both nonlinear and nonlocal effects associated with a 
dj.2_y2-wave order parameter [2^ and does not strictly 
reflect the field dependence of the superfluid density. In 
sample TWl which constitutes the most complete data 
set for YBa2Cu3 6.95, r-o(O) = 120.7 A, ^ = 0.075 T"! 
and 7 = 1.82 T^^. On the other hand, in the twinned 
YBa2Cu3O6.60 sample of Ref. ||l^ we find that ro(0) = 
156.3 A, (3 = 0.094 T^^ and 7 = 1.77 T^^. The solid 
curves in Fig. 4 represent Eq. (|2|) with the corresponding 
values of ro(0), /? and 7. These extrapolations suggest 
that the shrinking of the vortex cores saturates at large 
H. Note that Eq. (|^) can be written strictly in terms of 
the intervortex spacing L (since L oc Vi?) . 

Maggio-Aprile et al. attributed the two peaks ob- 
served in the spectrum for tunneling into a vortex core 
of YBCO at T = 4.2 K and iJ = 6 T, to the lowest 
bound quasiparticle energy level -E1/2 = 5.5 meV. Using 
the formula Ei/2 — 2^Aq/ Ep from Ref. |Q] and taking S^ab 
at low T to be equivalent to the BCS coherence length 
io = hvf/TTAo, gives U = {2h^ /men^Ei/^y^^ w 17 A. 
This result (« ro) is plotted in Fig. 4. The agreement 
with the extrapolated curve from our ^SR measurements 
|p5| is striking and raises the possibility that the vortex 
cores in YBa2Cu3 06.95 are conventional-like {i.e. they 
contain localized quasiparticle states). Several theoreti- 
cal studies {e.g. Q) have suggested that bound states 
may arise from a second component in the order param- 
eter, induced by spatial variations in the dj.2_j,2-wave 
component in the vicinity of a vortex core. If this is 
the case, the core expansion will be accompanied by the 
formation of numerous bound states, which should be de- 
tectable by STM. Thus far, STM has not been used to 
probe the vortex structure at low magnetic fields in the 
HTSs. 

In conclusion, we have observed a large increase in the 
size of the vortex cores in YBa2Cu3 6.95 at low magnetic 
fields, similar to the behaviour reported in the conven- 
tional superconductor NbSe2. The expansion of the cores 
appears to be a general property of superconductors in 
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the vortex state. The agreement with STM measure- 
ments based on a conventional treatment of the vortex 
core in YBa2Cu3 06.95, supports the existence of localized 
states in the core. 

We thank K. Machida for helpful discussions. This 
work was supported by NSERC and DOE grant DE- 
FG05-88ER45353. 



FIGURE CAPTIONS 

Figure 1. (a) The muon precession signal in 
YBa2Cu3 06.95 after field cooling to T = 5.8 K in a mag- 
netic field iJ = 0.498 T. The solid line is a fit described in 
the text, (b) Difference between the measured precession 
signal and the theoretical polarization function [in (a)] 
assuming S^ab = 54 A (solid circles) and S^ab = 20 A (open 
squares). For visual clarity, we have doubled the bin size 
used in (a) and not shown the error bars (~ ±0.003). 

Figure 2. (a) The ratio x^/NDF as a function of (a) 
^ab and (b) k (from the same fits) for H — 0.498 T (stars, 
NDF=1148) and H = 1.952 T (circles, NDF=1196) at 
r = 5.8 K. 

Figure 3. (a) The temperature dependence of tq at 
0.5 T (open circles) and 1.5 T (solid circles) in twinned 
YBaaCusOe.gs (TWl). (b) The magnetic field depen- 
dence of ro extrapolated to T — Q. The twinned crystals 
TWl and TW2 are shown as open circles and squares, 
respectively, whereas the detwinned crystal is denoted by 
solid triangles. 

Figure 4. The magnetic field dependence of vq in 
YBa2Cu3 06.95 (solid circles) and YBa2Cu3O6.60 (open 
triangles) extrapolated to T — 0. The solid square at 
H — 6 T is the value of ^ab ~ ''o deduced from the STM 
experiment |^ on twinned YBa2Cu307_5 at r = 4.2 K. 
The solid curves are explained in the text. 
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